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Objectives 


♦ Determine the suitability of various inviscid flux schemes for: 

• Acoustics problems 

• Directly resolving turbulence 

• Large-eddy simulations of turbulence 

♦ Provide guidance for the long-term evolution of our Loci-based 
CFD codes 

• Focus is on the compressible Navier-Stokes equations due to wide 
range of applicability 

• Emphasis in this paper is on the basic resolving characteristics of 
spatial discretizations and inviscid fluxes, not on LES sub-grid scale 
(SGS) models 

• Grid resolutions considered here are from DNS levels to fine-grid 
LES levels 



DNS 


LES/ILES 


HRANS-LES 


Turbulence Modeling Strategies 


Directly resolve all scales of turbulent motion. There is no 
eddy viscosity, and the intrinsic dissipation of the inviscid 
flux must be « laminar viscosity. 


Directly resolve largest scales of turbulent motion, and 
model the smaller scales using one of two approaches: 

A) classic LES - a subgrid-scale (SGS) turbulence model, 
requiring that the intrinsic dissipation of the inviscid flux is 
« SGS eddy viscosity, or B) ILES - where the intrinsic 
dissipation of a 2 nd order accurate inviscid flux 
approximately mimics the SGS eddy viscosity. 


Attempts to directly resolve turbulence only in regions 
with adequate grid resolution, otherwise turbulence is 
modeled. In LES regions the algorithmic requirements 
would be consistent with the above description. 


RANS 


Turbulence is entirely modeled. 


Governing Equations of Fluid Mechanics 


Navier-Stokes Equations in Terms of Conservative Variables: 


dU d(F, - F v ) d(G, - G v ) d(H,-H v ) 

H" H" H" 

dt dx dy dz 


p 


pu 


0 

pu 


puu + P 


Txx 

pv 

^ F i= - 

puv 

- F v = - 

"*"xy 

pw 


puw 


"*"xz 

pEtot 


puH tot 


UT XX + VT xy + WT XZ - q x 


Conservative Differencing Forms the Flux at a Half-Point: 
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Central Difference Schemes 
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CD-2 F| +1/2 - F avg (i,i+1 ) 

CD-4 F i+1/2 = (4/3)F avg (i,i+1) - (1/6)[F avg (i-1,i+1) + F avg (i,i+2)] 

CD-6 F i+1/2 = (3/2)F avg (i,i+1 ) - (3/1 0)[F avg (i-1 ,i+1 ) + F avg (i,i+2)] 

+(1/30)[F avg (i-2,i+1) + F avg (i-1 ,i+2) + F avg (i,i+3)] 

CD-8 F i+1/2 = (1 6/1 0)F avg (i,i+1 ) - (4/1 0)[F avg (i-1 ,i+1) + F avg (i,i+2)] 

+(8/1 05)[F avg (i-2,i+1 ) + F avg (i-1 ,i+2) + F avg (i,i+3)] 
-(1/140)[F avg (i-3,i+1) + F avg (i-2,i+2) + F avg (i-1 ,i+3) + F avg (i,i+4)] 
CF-2 F i+1/2 = (3/2)F avg (i,i+1) - (1/4)[F avg (i-1,i+1) + F avg (i,i+2)] 


F avg (ii,i2) = y 2 ( Pi1 + p i2 ) y 2 ( Ui1 + u i2 ) y 2 (v M + v i2 ) + y 2 (P M + p i2 ) 

where V = [1 , u, v, w, H tot ] and P = [0, P, 0, 0, 0] 

As shown by Pirozzoli (2010), this F avg reproduces the skew-symmetric scheme of 
Kennedy and Gruber (2008). Note if F avg = (F i1 + F i2 ), then standard divergence form 
central differencing formulas result. 


Upwind Schemes 
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Roe Flux using “left” and “right” interpolations for the half-node 

Fj+i/2 = y^L + F i+1/2 , R ) - A roe | ( U i+ y 2jR - Uj +1/2jL ) 

where F j+1/2L and U i+y2L are formed from [p, u, v, w, P] i+1/2 L 
and F i+1/2jR and U i+ y 2jR are formed from [p, u, v, w, P] i+1/2jR 

Primitive variable reconstruction for generic flow variable q> (no slope limiting) 


UB-1 cp i+1/2 , L = cp| 

UB-3 cp i+1/2iL = (-cp M + 5cp| + 2(p i+1 )/6 

UB-5 cp i+1/2 L = (4cpj. 2 - 26q> M + 94tPj + 54cp i+1 - 6<p i+2 )/1 20 

UB-7 cp i+1/2 L = (-6<Pi_ 3 + 50cpj_ 2 - 202cp M + 638cp| + 428cp i+1 - 76<p i+2 + 8cp i+3 )/840 

UF-2 cp i+1/2L = (-<p M + 4cpi + <p i+1 )/4 

CU-5 3(pi_, /2L + 6cp i+ i /2iL + cp; +3/ 2 L = ((Pm + 19(p, + 10cp i+1 )/3 


The cp i+ i /2jR values are formed from a flipped interpolation 
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Modified Scaled Wavenumber, k r ' 


Fourier Analysis Characteristics 


Real Part of Modified Scaled Wave- 
number: Dispersion (Phase) 



Imaginary Part of Modified Scaled 
Wave-number: Dissipation 



Scaled Wavenumber, k 


♦ The dispersion (phase) characteristics of each UB scheme are the same as the CD scheme of 
one higher order of accuracy 

♦ UB-1 and CD-2 have a limited range of phase accuracy, while CU-5 has a relatively wide range 

♦ Central schemes are non-dissipative, while all upwind schemes experience increasing 
dissipation error at higher wavenumbers 

♦ Higher-order upwind schemes have smaller dissipation error at lower wavenumbers than lower- 
order upwind schemes 
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1-D Acoustic Standing Wave Problem 


Air at P 0 = 101325 Pa, T 0 = 25 °C, R = 287 J/Kg-K, p 0 = 1.18413 kg/m 3 , 

a 0 = 346.117 m/s 


u(x,t=0) = U 0 cos(nx/l), where 0 < x < 2ttI, I = 1 m, U 0 = 0.1 m/s, 
and n is the desired initial number of wavelengths 


Grid: 128 grid points, uniform spacing, periodic B.C. 


Assuming isentropic flow and small disturbances, the Euler equations can be 
simplified to the linearized equations of gas dynamics. The resulting exact 

solution is 

U(x,t) = 0.5U 0 (cos(n(x - at)/l) + cos(n(x + at)/l)) = U 0 cos(nx/l) cos(nat/l) 

P(x,t) - P 0 = 0.5p 0 a 0 U 0 (cos(n(x - at)/l) - cos(n(x + at)/l)) = 

p 0 a 0 U 0 sin(nx/l) sin(nat/l) 


4-stage Runge-Kutta time advancement, CFL = 1.0, 0.5 or 0.25 
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Pressure History at First Pressure Antinode 

n = 8, 16 PPW, k = tt/8 


Central Difference Schemes Upwind Schemes 



♦ Some dispersion error noticeable for CD-2 scheme 

♦ Dissipation is already severe for UB-1, barely noticeable for UF-2 and UB-3 



Pressure History at First Pressure Antinode 

n = 16, 8 PPW, k = tt/4 


Central Difference Schemes Upwind Schemes 



♦ Dispersion error for CD-2 scheme now more significant 

♦ Some dispersion error for CF-2 noticeable in opposite direction 

♦ Dissipation even worse for UB-1, and now significant for UF-2 and UB-3 


10 


Pressure History at First Pressure Antinode 

n = 32, 4 PPW, k = tt/2 


Central Difference Schemes Upwind Schemes 



♦ Dispersion error severe for CD-2, significant for CD-4, and noticeable for 
the other central schemes 

♦ Dissipation severe for UB-1, UB-3 and UF-2, and significant for all other 
upwind schemes 

♦ The least dissipative is CU-5, which also has the best dispersion accuracy 
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Taylor-Green Vortex Problem 


Air at P 0 = 7271 Pa, T 0 = 25 C C, R = 287 J/Kg-K, p 0 = 0.0849723 Kg/m 3 , 
a 0 = 346.117 m/s, U 0 = 0.1 aO = 34.6117 m/s, |j 0 = Sutherland relation, 

Re = p 0 U 0 l/Mo = 1600 


u(x,y,z,t=0) = U 0 sin(x/I) cos(y/l) cos(z/l), 
v(x,y,x,t=0) = -U 0 cos(x/l) sin(y/l) cos(z/l), 

w(x,y,x,t=0)=0, 

p(x,y,z,t=0) = p 0 + (p 0 U 0 2 /16) [cos(2x/l) + cos(2y/l)][cos(2z/l) + 2] 

p(x,y,z,t=0) = p/RT 0 

where -ttI < x < ttI, -ttI < y < ttI, -ttI < z < ttI, I = 0.01 m 

Grids: 256 * 256 x 256, 192 x 192 x 192 and 128 x 128 x 128 grid points, 

uniform spacing, periodic B.Cs. 


4-stage Runge-Kutta time advancement, CFL » 0.45 on 256 x 256 x 256 grid 


DNS results have been calculated by Brachet et al. (1983) and van Rees et al. 
(2011). Also see http://www.public.iastate.edu/zjw/hiocfd.html. 
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Iso-surfaces of Q-criterion Q = 0.1 (U 0 /l) 2 
Taylor-Green Vortex Problem at Re = 1600 



t* = 12 t* = 1 6 t* = 20 

♦ Simulation run using CD-8 scheme on 256 x 256 x 256 point grid 

♦ Iso-surfaces colored by velocity magnitude U/U 0 

♦ t* is a normalized timescale tU 0 /l 
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Mean Kinetic Energy, E k,m /V ; 


Mean Kinetic Energy History 
256 x 256 x 256 Point Grid 


Central Difference Schemes 


Upwind Schemes 




♦ The central schemes generally follow the DNS kinetic energy history very 
well, though CD-2 differs slightly toward the end 

♦ The UB-1 scheme rapidly dissipates all kinetic energy in the flow 

♦ UB-3 and UF-2 are much closer to the DNS, though some differences are 
evident. 


♦ The remaining upwind schemes follow the DNS very closely 


14 


Mean Kinetic Energy Dissipation Rate, s l/v; 



Measured Mean Kinetic Energy Dissipation 
Rate History, 256 x 256 x 256 Point Grid 


Central Difference Schemes Upwind Schemes 

CO O CO o 



♦ The central schemes generally follow the DNS kinetic energy dissipation 
history very well, though CD-2 misses the peak at t* = 9 

♦ The UB-1 scheme has such a high dissipation rate that it is not shown 

♦ UB-3 and UF-2 are much closer to the DNS, though both are overly 
dissipative through t* = 8, and peak earlier at a lower value 

♦ The remaining upwind schemes follow the DNS very closely 
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Mean Enstrophy, Q m I 2 / V 2 


Mean Enstrophy History 
256 x 256 x 256 Point Grid 


Central Difference Schemes 


Upwind Schemes 




♦ The central schemes generally follow the DNS enstrophy history very well, 
though CD-2 misses the peak at t* = 9 

♦ The UB-1 scheme rapidly dissipates all enstrophy in the flow 

♦ The remaining upwind schemes also all fail to match the DNS, with UF-2 
having the worst agreement, and CU-5 the best 
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Iso-contours of Vorticity Magnitude Ml/U 0 
256 X 256 X 256 Point Grid at t* = 8 



Note: DNS was run 
on a512 x 512 x 512 
point grid 


DNS CD-2 


CD-4 


CD-6 





Contour Levels: 
1,5, 10, 20 and 30 


CD-8 CF-2 UB-1 UB-3 



x = 0 plane 
ttI/ 2 < y < ttI 
ttI/2 < z < ttI 


UB-5 UB-7 


UF-2 


CU-5 


♦ Dissipation severe for UB-1, very noticeable for UF-2 and UB-3 


♦ Good results for other schemes, though lack of sharpness for CD-2 
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Mean Kinetic Energy, E km / V; 


Mean Kinetic Energy History 
128x 128x 128 Point Grid 


Central Difference Schemes 


Upwind Schemes 




♦ The central schemes generally follow the DNS kinetic energy history 
reasonably well, though CD-2 differs slightly toward the end 

♦ The UB-1 scheme rapidly dissipates all kinetic energy in the flow 

♦ UB-3 and UF-2 noticeably diverge from the DNS 

♦ UB-5, UB-7 and CU-5 follow the DNS reasonably well 
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Mean Kinetic Energy Dissipation Rate, s I / V; 


Measured Mean Kinetic Energy Dissipation 
Rate History, 128 x 128 x 128 Point Grid 


Central Difference Schemes Upwind Schemes 

CO O CO o 



♦ The central schemes generally follow the DNS kinetic energy dissipation history 
reasonably well, though CD-2 peaks higher and earlier than the others 

♦ The UB-1 scheme has such a high dissipation rate that it is not shown 

♦ UB-3 and UF-2 have a significantly higher dissipation rate than the DNS for t* = 0-8, 
and have a lower peak earlier in time 

♦ UB-5, UB-7 and CU-5 do much better, but none capture the peak of the DNS 
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Mean Enstrophy, Q m I 2 / V 2 


Mean Enstrophy History 
128x 128x 128 Point Grid 


Central Difference Schemes Upwind Schemes 



♦ The central schemes generally follow the DNS enstrophy history very well, 
though CD-2 misses the peak at t* = 9 

♦ The UB-1 scheme rapidly dissipates all enstrophy in the flow 

♦ The remaining upwind schemes also all fail to match the DNS, with UF-2 
having the worst agreement, and CU-5 the best 


20 



Turbulent Channel Flow Problem 


Air at P 0 = 25331 Pa, T 0 = 25 °C, R = 287 J/Kg-K, p 0 = 0.296033 Kg/m 3 , 
a 0 = 346.117 m/s, u m = 44.44 m/s, |j 0 = Sutherland relation, 

Re 2 h = Po^m^h/Mo = 14300 


u(x,y,z,t=0) = u c0 [1 - (y/h) 8 ] + U 0 tt cos(2x/h) sin(TTy/h) sin(2z/h), 
v(x,y,x,t=0) = -U 0 sin(2x/h) [1 + cos(TTy/h)] sin(2z/h), 
w(x,y,x,t=0)= -(U 0 tt/ 2) sin(2x/h) sin(TTy/h) cos(2z/h) 
where 0 < x < 2ttIi, -h<y<h, 0<z< nh, h = 0.01 m, u c0 = 50 m/s, U 0 = 5 m/s 


Grids: 128 * 129 x 128, 96 x 129 x 96, 64 x 129 x 64 and 48 x 129 x 48 grid 
points, uniform spacing and periodic B.Cs. in x- and z-directions, 
nonuniform spacing and viscous wall B.C. in y-direction 


4-stage Runge-Kutta time advancement, CFL » 2.45 based on wall-normal 

grid spacing at wall and acoustic wavespeed 


Survey of experimental results reported by Dean (1978), and extensive DNS 
studies [Kim et al. (1987) and Moser et al. (1999)] have also been performed 
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Iso-surfaces of Q-criterion Q = 0.1 (u m /h) 2 
Turbulent Channel Flow Problem at Re 2h = 14300 




t* = 30 


♦ Simulation run using CD-8 scheme on 128 x 129 x 128 point grid 

♦ Iso-surfaces colored by velocity magnitude U/u m 

♦ t* is a normalized timescale tu m /2TTh 
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Streamwise Velocity, 


Streamwise Velocity Profiles 
128 x 129 x 128 Point Grid 


ID 


♦ Time- and space-averaged results are normalized by wall-variables 

♦ All of the central schemes exhibit excellent agreement with the DNS, and merge 
smoothly from the laminar sublayer to the log layer 

♦ The UB-1 scheme (not shown) completely laminarizes the flow 

♦ UB-7 and CU-5 also exhibit very good agreement with the DNS 

♦ UB-5, UB-3 and UF-2 are progressively worse, and increasingly rise above the log layer 




10 ° 10 1 _ + 10 2 
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Central Difference Schemes 


Laminar Sublayer: u + = y + / 

Log Layer: u + = 2.44*ln(y + ) + 5.2 , 

_ DNS, Re =395 , 
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Upwind Schemes 


m — — | • I ' — ! 1 1 1 1 'I 

Laminar Sublayer: u + = y + / 

Log Layer: u + = 2.44*ln(y + ) + 5.2 / 


RMS Velocities 


RMS Velocity Profiles 
128 x 129 x 128 Point Grid 



Central Difference Schemes 


0.2 0.4 0.6 0.8 1 

Distance Across Channel Half-Width, |y - yj/h 
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Upwind Schemes 


0.2 0.4 0.6 0.8 

Distance Across Channel Half-Wdth, |y - yj/h 
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8 30 
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CD 

« 20 


♦ Time- and space-averaged results are normalized by wall-variables 

♦ The central schemes also show generally excellent agreement with the DNS RMS 
velocities, though CD-2 is not quite as good for u’ + 

♦ The UB-1 scheme (not shown) completely laminarizes the flow 

♦ UB-7 and CU-5 also exhibit very good agreement with the DNS 

♦ UB-5, UB-3 and UF-2 are progressively worse, with u ,+ peaking at higher values 
further away from the wall, while the peaks of v ,+ and w ,+ decrease 
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Turbulent Shear Stress, -u'v'/u 


Turbulent Shear Stress, -u’v’/u , 2 
128 x 129 x 128 Point Grid 


Central Difference Schemes Upwind Schemes 



♦ Time- and space-averaged results are normalized by wall-variables 

♦ All of the central schemes exhibit excellent agreement with the DNS 

♦ The UB-1 scheme (not shown) completely laminarizes the flow 

♦ UB-7, CU-5 and even UB-5 also exhibit very good agreement with the DNS 

♦ UB-3 and UF-2 are progressively worse 
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Streamwise Velocity, 


Streamwise Velocity Profiles 
48 x 129 x 48 Point Grid 


Central Difference Schemes Upwind Schemes 



♦ Time- and space-averaged results are normalized by wall-variables 

♦ CD-6, CD-8 and CF-2 are consistently offset below the log layer 

♦ CD-4 and CD-2 are closer to the log layer but now have the wrong slope 

♦ UB-1, UB-3 and UF-2 (not shown) completely laminarize the flow 

♦ CU-5, UB-7 and UB-5 are progressively offset above the log layer 


RMS Velocities 


RMS Velocity Profiles 
48 x 129 x 48 Point Grid 


Central Difference Schemes 


Upwind Schemes 



-n | i r 

DNS, Re =395 

CD-2 

CD-4 

CD-6 


0.2 0.4 0.6 0.8 

Distance Across Channel Half-Width, |y - yj/h 



♦ Time- and space-averaged results are normalized by wall-variables 

♦ None of the central difference schemes match the DNS RMS velocity profiles very 
well, with CD-2 having the worst agreement 

♦ UB-1, UB-3 and UF-2 (not shown) completely laminarize the flow 

♦ CU-5, UB-7 and UB-5 are progressively worse in predicting the DNS results, with u ,+ 
peaking at higher values further away from the wall, while the peaks of v ,+ and w ,+ 
decrease 
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Turbulent Shear Stress, -u'v'/u 


Turbulent Shear Stress, -u’v’/u, 2 
48 x 129 x 48 Point Grid 


Central Difference Schemes Upwind Schemes 



♦ Time- and space-averaged results are normalized by wall-variables 

♦ All of the central schemes exhibit good agreement with the DNS, though 
some differences are more evident for CD-2 than the other schemes 


♦ UB-1, UB-3 and UF-2 (not shown) completely laminarize the flow 

♦ CU-5, UB-7 and UB-5 are in progressively worse agreement with the DNS 
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Normalized Skin Friction Coefficient, c f /c 



Effect of Spanwise Grid Spacing on Average 
Skin Friction Coefficient, c f /c f dean 


Central Difference Schemes Upwind Schemes 



♦ C f dean = 0.073 Re 2h ' 0 25 = 6.68 x 10 -3 for these cases 

♦ All central schemes provide a good prediction of c f for Az + <10, while at 
larger values of Az + the higher-order schemes overpredict c f . CD-2 
underpredicts c f on the coarsest grids. 

♦ CU-5 and UB-7 yield an excellent prediction of c f for Az + <10, and a 
reasonable prediction (within 10%) for for Az + <13. UB-5 may be good for 
for Az + » 5. UB-3 and UF-2 are very poor in predictive accuracy. 
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Summary 


♦ Five different central difference schemes, based on a conservative 
differencing form of the Kennedy and Gruber skew-symmetric scheme, 
were compared with six different upwind schemes based on primitive 
variable reconstruction and the Roe flux 

♦ These eleven schemes were tested on a one-dimensional acoustic 
standing wave problem, the Taylor-Green vortex problem and a turbulent 
channel flow problem 

♦ The central schemes were generally very accurate and stable, provided the 
grid stretching rate was kept below 10%. As near-DNS grid resolutions, the 
results were comparable to reference DNS calculations. At coarser grid 
resolutions, the need for an LES SGS model became apparent. There was 
a noticeable improvement moving from CD-2 to CD-4, and higher-order 
schemes appear to yield clear benefits on coarser grids 

♦ The UB-7 and CU-5 upwind schemes also performed very well at near-DNS 
grid resolutions. The UB-5 upwind scheme does not do as well, but does 
appear to be suitable for well-resolved DNS. The UF-2 and UB-3 upwind 
schemes, which have significant dissipation over a wide spectral range, 
appear to be poorly suited for DNS or LES. 
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